Scaling of Nonlinear Longitudinal and Hall Resistivities near the Vortex Glass 

Transition 
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We show that the longitudinal current-voltage characteristics of superconductors in mixed state 
have the general form of extended power law. Isotherms simulated from this nonlinear equation fit 
the experimental I-V data of Strachan et al. [ Phys. Rev. Lett. 87, 067007 (2001)]. We determine 
the average pinning force in the flux creep and strong pinning regime and discuss both the puzzling 
scaling behavior pxy oc j. and a recently found new scaling relationship of nonlinear Hall resistivity 

pxy{T). 
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The vortex motion in the superfluid electrons has 
presented a persistent problem in the physics of high- 
temperature superconductors(HTS). One of the most 
puzzling phenomena is the sign change that has been 
observed in the Hall effect in the superconducting state 
in most HTS and some conventional superconductors . 
Another highly intriguing issue is the power law scal- 
ing relation between the Hall resistivity pxy and longitu- 
dinal resistivity p^x Qi HI • Experiments show that the 
Hall effect at least depends on two factors: the doping 
level and the vortex pinning 0, 0, 01 • Based on the 
hypothesis of a vortex glass transition 0, Dorsey and 
Fisher Q reasonably explained the puzzling scaling be- 
havior Pxy oc [Pxx{T)f with /5 = 1.7 ± 0.2 observed by 
Luo et al. 0. The Hall effect itself was attributed to 
a " particle- hole" asymmetry and the exponent (3 was re- 
lated to a specially chosen particle-hole asymmetry ex- 
ponent A ~ 3. They predicted "The nonlinear Hall 
electric field Ey should exhibit universal scaling and, 
right at the transition, should vanish with a universal 
power of the current density J^:". This prediction has 
been verified by Woltigens et al. 0] with an experimen- 
tally determined Hall-related exponent A = 3.4 ± 0.3. 
Furthermore, they found wide range scaling behavior 
Pxy — A[pxx{Jx,T)]'^'^'^^'^ in consistency with the model 
proposed by Vinokur, Geshkenbein, Feigel'man and Blat- 
ter(VGFB) j^]. However, recently the widely accepted 
experimental evidence for vortex glass theory has 
met serious doubt ^Jl- Strachan et al. show wide 
range accurate isothermal I-V measurements over 5 or 
6 decades on a high quality 2200 A thick YBa2Cu307_a 
film with Tc « 91.5K and transition width about 0.5K 
PH and find although the I-V isotherms measured in 
a magnetic field can be collapsed into scaling functions 
proposed by Fisher et al. Q as is widely reported in 
the literature, these excellent data collapse can also be 
achieved for a wide range of exponents and glass tem- 
peratures Tg as demonstrated in their Fig. 2(a), 2(b) and 
2(c) _llj. Since the critical temperature Tg cannot be 
determined uniquely, the correctness of the vortex glass 
picture has to be reinvestigated. In this work, we show 
that the longitudinal isothermal current-voltage charac- 
teristics have the general form of extended power law |12| . 



This remarkable feature naturally explains the insightful 
new findings of Strachan et al. jllj] as well as the recently 
observed scaling relationship of pxy{T). 

Thermally activated flux motion can be considered as 
the sequence of thermally activated jumps of the vor- 
tex segments or vortex bundles between the metastable 
states generated by disorder. Every elementary jump is 
viewed as the nucleation of a vortex loop, and the mean 
velocity of the vortex system is determined by the nucle- 
ation rate 0, [l3| 



V (X exp {-SF/kT) . 



(1) 



Here SF is the free energy for the formation of the crit- 
ical size loop or nucleus which can be found by means of 
the standard variational procedure from the free energy 
functional due to the in-plane displacement u (z) of the 
moving vortex during loop formation 
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where fL = <&oJ x Ez/c is the Lorentz force due to ap- 
plied current J and is the viscous drag force on vortex, 
= — ?7Vvortex, with Vvortex = du/dt and viscous drag 
coefficient « ($0-802) / (pnC^), Pf is the flux flow resis- 
tance of a pinning free mixed state w , as derived 



by Bardeen and Stephen 14]. 

Our equation jSJ differs from the corresponding one 
used in previous vortex glass theories on that we con- 
sidered both the Lorentz force Tl and the viscous drag 
force iri on the vortex loop during its formation while in 
present vortex glass theories ordy bare Lorentz force fL 
has been taken into account 0, 0, 0] • 

As shown in 16], the barrier energy found by the con- 
ventional vortex glass model 7] has the general form with 



2 



that of collective pinning model and Boson glass model 



SF^Fioop[L''iJ)]=U{J)^Uo{j-r, (4) 

with L*{J) the critical size of vortex loop nuclcation. 
Since the viscous drag force f,, has been taken into ac- 
count in Eq. Q , now instead of Eq. Q one must have the 
barrier energy 



SF = U{Jp) « C/o(^)^ 



(5) 



which implies a current- voltage characteristic of the form 



Eij) = pfjcxp[-^{^r] 



(6) 



where Uq is a temperature and field dependent charac- 
teristic pinning energy related to the stiffness coefficient 
and Jo is a characteristic current density related to Uq, 
fi is an numerical exponent. 

Considering the real sample size effect we note that the 
barrier energy for vortex loop excitation cannot exceed 
the free energy value needed for generating a loop with 
real sample size D defined as 

SF = UiJ) < U{Jd) = Fioop[L*iJ = Jd) = D]. (7) 

Therefore, the free energy barrier for vortex loop exci- 
tation never diverges even at vanishing applied current J 
as predicted by the conventional vortex glass theory • 

Replacing Jp in Eq.® with the expression Eq.©, and 
taking the logarithm of its both sides, one obtains 
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,(8) 



where J^f = E{JD)/pf is much smaller than Jjj and 
we define h ^ -['"(■/W./../);;"-in(J/./,)^/^] ^ 

proximation (l-t-/i)~^ w 1 — h for \h\ < 1, one finally 
obtains a general normalized form of the current-voltage 
characteristic as 



2/ = a; exp [-7 (1-f ?/ - xf] 



(9) 



where 



kT \Jd 
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2Jd' 
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EjJ) 

2pfJD 



, p = /i. 



Equation® can also be nearly equivalently expressed 
in an extended power law form 



y 



(10) 



with n is a parameter determined by the condition that 
Eqs.Q and meet at the inflection points of Iny ~ 
Inx curves for Eq.®. 

In Fig. 1(a) we show the numerical solutions of Eq.® 
and Eg. 1)10(1 for comparison. 

The complex longitudinal experimental current- 
voltage data are usually described phenomenologically 
with a form [13 



EiJ) = Jp^^iT,B,J) = Jpfe 



-U^ff(T,B,.J)/kT 



(11) 



Equation® implies that effective barrier can be explic- 
itly expressed as 

Ueff{T,B,J) = U,{B,T)F[J/UB,T)l (12) 

with Jc{B,T) = 2Jd{B,T). 

Incorporating it into the commonly observed scaling 
behavior of magnetic hysteresis M{H) in superconduc- 
tors, it can be shown that Uc{B,T) and Jc{B,T) in 
Ea. ((12|l must take the following forms 

Uc{B,T) = ■^{T)B'^ oc [r*(B) -r]*B", 
J^(B, T) = X{T)B"' cx [T*{B) - T]°'B"', (13) 

with T*{B) being the irreversibility temperature where 
tilt modules and stiffness of vortex matter vanish. Con- 
sidering pf oc T, one finds from Ea. ((12|l the follow- 
ing temperature dependent relations for 7 and I — V 
isotherms of optimally doped YBCO samples in a given 
magnetic field as 



7(T) = 7o(r*-r)Vfcr, 

/ oc xJd{T) oc x(T* -T)", 

V o^yJD{T)pf o,y{T* -T^T. 



(14) 



The combination of equations ® and 114() enables 
us to reproduce the experimentally measured I — V 
isotherms of Strachan et al. ini. In Fig. 1(b) we show the 
comparison, where the parameters used are a = 3.0, 5 = 
0.5, p — 1.75,70 = 3.21 x 10~^^ and the irreversibility 
temperature T*{B = AT) is assumed SSJC similar to the 
Fig.l in Ref.ini- 

On the basis of nearly 100 simulated I ~V isotherms 
from to 88X with temperature intervals of 0.2Jir, we 
obtain three different scaling data collapse analysis with 
assumed critical temperature Tg of %1K JbK and 70K 
respectively. In Fig. 2 we show the comparison of our 
result with the data collapses of measured I — V curves 
of Strachan et al. in Ref. and find fair agreement. 

The analytical extended power law equations ®, 110|) 
and (|12|) enable one to resolve the nontrivial problem of 
determining the positive scaling function T{v) between 
the average pinning force (Fp)t and vortex velocity vl 
used in the VGFB model as well as by Wang, Dong 
and Ting (WDT) (iTf in the form (Fp)t = -r(t;L)vL. 
Thus, the analytical equation for the nonlinear Hall re- 
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FIG. 1: (a)Comparison of solutions of Eq.@(solid symbols) 
and Eg. mUll fcurves). 

(b) Comparison of our simulated I-V isotlierms(lines) with the 
experimental data of Strachan et al. flU (solid symbols). 



sistivity can be derived as 

2 

A*„,p„i?{(l + C)%^ -2C^}, 

Pf 



(15) 



with the longitudinal nonlinear resistivity defined in 
Eas. lfTT|l and ((T^ . p^x — Pfe^^"!^^ ■ The parameter 
C is defined as C = C(l — H jHc-i) with H the aver- 
age magnetic field over the vortex core and C is describ- 
ing the contact force on the surface of vortex core [l7| . 
/3o = fJ'mHc2 with = Te/m the mobility of the charge 
carrier and Hc2 = $o/27r^^ the upper critical magnetic 
field . 

This equation is based on the conventional theories of 
Bardeen and Stephen as well as Nozieres and Vinen 
lis"! for pinning free motion of vortex which predict that 
the Hall effect stems from quasinormal core and hence 
has the same sign as in the normal state. However, the 
microscopic theories of the vortex dynamics in clean su- 
perconductors ^3 predict that even in the absence of 
pinning, the vortex contribution to Hall current may have 
the sign different from the sign of Hall effect in the nor- 
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FIG. 2: •: Collapse of the simulated I-V isotherms repro- 
duced from Eq.lO and Eq.l^lJ using various critical param- 
eters compared with 

o: the experimental I-V isotherms for a 2200A YBCO film 
in 4T measured by Strachan et al. flU . The conventional 
analysis is shown in (a). 



mal state. In latter case, increasing pinning can remove 
the Hall anomaly and even result in a second sign rever- 
sal of An/ observed in some strongly anisotropic materials 

With consideration of this factor. Hall resistivity equa- 
tion (|15fl should be modified to the form 

P.y{T)^PrnPnB-^^{il + C)^ - 2C^}, (16) 
An(-L) Pf Pf 

where As (T) and A„ (T) are the Hall force coefficients of 
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FIG- 3: Scaling of the anomalous pxy 



T curves of 



different kinds of cuprate superconductors in the form of 
logio\p:cy/ pm \ ~ lopio((T - Tq) / {Tm " Tq)) . Therein (a) rep- 
resents YBaCusOr 1,^0]; (b) represents Bi2Si2CaCu2 08-)-a; 
[2H : (c) represents HgBa2CaCu2 06-f 5 (J]; (d) represents 
Tl2Ba2CaCu2 08 H- 



the sample in pinning free superconducting and normal 
states respectively ^19]. 

Besides the well-known puzzling scaling behavior 
Pxtj ~ Pxx' recently we find a new universal scaling rela- 
tion of the form 



log 



10 



Pxy 
Pm 



T-Tq 

Tm. — Tn 



(17) 



where pm is the first extreme value of Hall resistivity 
from low temperature side and Tm is the temperature 
corresponding to pm- Tq corresponds to the temperature 
where the Hall resistivity become first measurable at the 
low temperature side. 



This scaling relationship is consistent with the gen- 
eral Hall resistivity eauation Hl()|l . In the case of YBCO 
(Fig.3(a) H 113), the factor A,{T)/An{T) > 0, we 
see the Hall anomaly at temperatures around Tm, on 
the other hand, for the strongly anisotropic materials 
HBCC0(Fig.3(c) 11) and TBCC0(Fig.3(d) [H), we see 
the double sign reversal of pxy{T), perhaps owing to that 
As{T)/An{T) < 0, as mentioned by Kopnin and Vinokur 



In conclusion, we demonstrate that the longitudinal 
E{ J) response have the general form of extended power 
law. This nonlinear behavior naturally explains the scal- 
ing of longitudinal and Hall resistivities near the vortex 
glass transition. 
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